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ABSTRACT 


The relative distance between various Antisubmarine Warfare vehicles 
is an important parameter for the evaluation of sensor equipment. This 
thesis presents a computer model, utilizing a time dependent error navi- 
gation system, Omega, for the determination of this relative distance 
parameter. A methodology for predetermining a relative distance error 
distribution for a specific area by simulating various Antisubmarine 
Warfare vehicle tracks from known navigational errors has been estab- 
lished. Attempts to smooth estimated tracks of Antisubmarine Warfare 
vehicles which receive the same time dependent navigational information 
simultaneously was found not to be advisable as this disturbs the correla- 


tion of navigation error between the vehicles. 
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Section One 
Introduction to the Problem 
1-1 General Discussion of the ASW Exercise 

The fleet's AntiSubmarine Exercises are an excellent example of the 
operational environment in which test and evaluation of weapon systems can 
be conducted. These exercises are necessary, for the sophistication and 
cost of the present state-of-the-art of military weapon systems requires 
that their capability be determined in this test and evaluation phase. 

The ASW exercises not only provide the opportunity to observe the 
new system but also to train personnel and determine the effectiveness of 
the existing tactics and equipment. However, this method of evaluation 
is expensive which places additional emphasis on the essentiality of ob- 
taining the best results which will in fact reflect the true capability 
of the weapon system. 

1-2 Problem 

One of the primary purposes of an ASW exercise is to provide data to 
evaluate the potential of detecting devices of the ASW forces. The inter- 
play of the different vehicles involved (i.e. aircraft, surface ships and 
Submarines) in an ASW exercise create a unique problem in that not only 
are their exact positions never established but also confirmation of de- 
tection by actual sightings are not always possible. This situation 
breeds the continually arising question, 'Were any two vehicles within de- 
tection range of the detecting sensors?" 

This question is hopefully resolved during the formal reconstruction 
of the ASW exercise which takes place after its termination. The person- 
nel usually represented at these sessions are a mixture of those directly 


involved, along with the necessary analyst and administrative personnel. 


One of their first tasks is to determine the tracks of the various vehicles 
by reconstructing them from the respective navigators' logs. The present 
reconstruction procedure is spelled out in great length which results in 

a process that is tedious and time consuming. (1 | Unfortunately, even after 
these tracks have been reconstructed, serious doubt exist to their accu- 
racy. The cause of this doubt does not arise from the reconstruction pro- 
cedures but basically in its source of navigational information--the 
navigator's log. 

The navigator's log is composed of information gathered from available 
navigational equipment on board the respective vehicles. The navigational 
errors present in the log are the result of an accumulation of navigation- 
al equipment and human errors. The relative distance between any two ve- 
hicles in question is of primary interest in this reconstruction session. 
The error in this relative distance between any two vehicles (i.e. see 
figure 1) results from a combination of their separate inherent naviga- 
tional errors. Knowledge concerning this relative distance error is 
presently not well formulated. But this relative distance error is one 
of the parameters necessary in analyzing the effectiveness of the detec- 
tion equipment. 

1-3 Objective 

General Electric TEMPO, Santa Barbara, California is presently under 
contract with the United States Navy to investigate the problems intro- 
duced by the relative distance errors and to propose means of alleviating 
or reducing them. Smoothing techniques are being considered which would 
attempt to reduce the existing relative errors by smoothing the estimated 
positions found from present navigation information. The employment of a 


new type of navigational aid called OMEGA, which is in the operational- 
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evaluation stage, is under consideration as the primary source of naviga- 
tional information. 

Tt is therefore our objective to investigate the aforementioned alter- 
natives and recommend a procedure that will aid considerably in the track 
reconstruction problem. Our recommendations will focus on two specific 
aspects, first a faster and more effective method of reconstructing the 
various ASW vehicles tracks and second to provide a means by which one 
could conduct a statistical error analysis of the tracks. 

Hopefully our solution would effectively free the personnel at the 
reconstruction session from doing the actual track reconstructing them- 
selves, thus permitting them to devote more time and energy towards strict- 
ly analyzing the exercises. 

1-4 Approach 

The first step is to tackle the basic underlying problem, that of 
errors introduced by all the various methods in which the navigation in- 
formation is presently obtained. It cannot be overemphasized that if er- 
rors must exist, it is highly desirable that these errors be dependent in 
respect to the navigational systems of both vehicles, when relative accu- 
racy is of interest. It follows that one should then concentrate on find- 
ing the optimal combination of navigational equipment required to provide 
this relative distance information. Clearly the minimum navigational in- 
formation necessary is the position of each vehicle at a given time. How- 
ever, in order that analysis be conducted, a sequence of such points must 
exist so that a continual picuture of the interplay is displayed. The 
time interval between the sequence of points depends on the speed at which 
the vehicles are traveling. The time interval between aircraft may be 


one minute while for submarine and surface vessels every five minutes may 
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be sufficient depending on the detail desired of the reconstructed tracks. 
This sequence of points is then basically a "dead-reckoning" track between 
navigational fixes for each vehicle. 

The present methods employed to obtain the DR tracks are to begin at 
a fix obtained from a navigation source and apply course and speed infor- 
mation as one's guide until another fix is established. This requires 
continuous monitoring of compasses and speed indicators if one expects to 
obtain the required accuracy desired for detection equipment evaluation. 
This process has not proven desirable even after highly sophisticated 
equipment has been employed to automatically plot the tracks. One may 
suggest that more fixes be obtained so that the DR time can be shortened 
thus resulting in better accuracy. At present this suggestion is sound in 
theory but is not practical. The reason being that short time intervals 
(e.g. less than fifteen minutes) between fixes for a majority of present 
Naval crafts is not easily obtained. 

The idea of receiving continuous fixes is intriguing, so let us pur- 
sue it further. A position can be obtained in two ways, first by knowing 
the range and bearing from a reference point. The second method is to 
find two lines of positions which intersect. Since reference positions 
(i.e. land masses) are non-existant the majority of the time, the inter- 
section of two lines of position appears to be the more desirable tech- 
nique. 

With this in mind, if one could monitor on board the ASW vehicles, 
some sort of navigational input signals that would establish their posi- 
tions, then a continuous fix could be realized. This could be accomplished 
by putting on board the vehicles recorders which record these signals and 


the corresponding time. This procedure would also take the place of the 
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navigator's log - eliminating human errors. 

Hopefully, if errors are introduced using this technique, they will 
be dependent to this technique and will not affect appreciably the relative 
accuracy between any two vehicles. 

The optimal navigation system selected to be used in conjunction with 
this technique is one that is not only compatible to a continuous readout 
but also possessing errors which result from a source common to all ve- 
hicles in the immediate vicinity, e.g. variations in the velocity of wave 
propagation. We propose that the Omega navigation system now in the op- 
erational-evaluation stage, which possesses these desired characteristics, 
be selected. However, any navigational system with these characteristics 
may be used. 

In fact, we now have a navigation system, consisting of Omega navi- 


gational information, which can be monitored continuously and simultan- 


eously by all vehicles. 
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Section Two 
Omega 
2-1 Introduction 

Omega is a Very Low Frequency (VLF) navigation system operating in 
the internationally allocated navigation band in the electromagnetic spec- 
trum between 10 and 14 Kc/s. (See Table 1 for further system character- 
istics). [2] As in Loran each pair of stations provides the navigator with 
a family of hyperbolic lines-of-position. The 5000 to 6000 nautical mile 
baseline between stations gives a global coverage capability with a maxi- 
mum of only eight stations. 

2-2 Propagation Medium 

Most modern radio navigation systems are designed to obtain optimum 
results from propagation medium involved. From this follows the fact that 
most navigation systems are propagation limited. Omega is no exception. 
Capabilities and operational procedures stem in large part from the prop- 
erties of VLF propagation. It should be noted that the propagation prop- 
erties to be considered are solely those at the Omega frequencies. Much 
data has been acquired at the basic 10.2 Ke/s frequency so definite con- 
clusions can be made. Omega is a simple system with range and accuracy 
based on the features of VLF propagation, second order effects are well 
understood although large. Easy to use corrections are prepared for use 
with Omega. 

The most apparent feature of an Omega VLF phase difference measurement 
as observed at some remote stationary point is a diurnal variation. While 
day-night transitions result in large changes in position, they are nor- 
mally quite regular and predictable. If the diurnal variation can be pre- 


dicted it need only be subtracted from observed data to yield a navigation 
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Receiver Specification (Ground Monitor) 

Frequency: 10.2 KC fixed or continuously 
Gunable. fromes: 2 teowr1>.2 KC. 

Sensitivity: 0.75 db noise figure. 

Differential Amplitude: 80 db. 

Input Impedance: 50 ohms. 

Detection: Synchronous. 


Instrumental Accuracy: 0.2% phase-difference 
of two noise-free signals. 


Tracking: Zero velocity error servos. 


Signal Conditions: Track reliably in signal-to- 
noise ratio of 1/10. 


Commutator Synchronization: Semi-automatic operation 
in signal-to-noise ratios 
LO Ly Oe 
Reference Oscillator: Long-time stability 1 part in 
108 per day, short time stability 
> parts in 109 per 30 minutes. 
Size: Height 57 *5", Width 19 45", Depth 17 45". 
Weight: 486 pounds. 


Power: 115 volts, single phase, 60 cycles, 
12 amperes. 


Table 1 


Omega Receiver Specifications 
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variate which, to a good approximation, varies quite regularly with posi- 
Elon. 
2-3 Accuracy 

Factors effecting the accuracy (i.e. the magnitude of the navigation- 
al error) to be obtained from a practical navigation system include: re- 
peatability of observations such as day-to-day variations at a particular 
site, accuracy with which observed readings may be associated with real 
terrain features, and degradation experienced in actual operations due to 
shipboard installation or procedures. 

Although Omega is an absolute system, relative accuracy is important 
in many applications. Because the limitations on Omega accuracy is the 
phase meawdiey of a transmission through long distances, phase perturb- 
ations of a signal from a remote transmitter at any two receiving sites 
in the same neighborhood will be much the same. Relative accuracy is 
therefore significantly better than absolute accuracy. It is also indica- 
tive of relative positioning capability with a task force. As separation 
distance between receivers decreases, the relative accuracy increases. 

The additional difference between the experimental and operational 
systems must be taken into account. The experimental system used short- 
term synchronization, that 18, propagational fluctuations on the path be- 
tween two Omega transmitters which were fed into the phase radiated by one 
of the transmitters. In the operational system, the synchronization will 
be achieved through the use of very stable oscillators, using a time con- 
Stant of one or more days. Thus the phase will be radiated at a constant 
value and corrected only slowly to keep it at the assigned value. 

The effect of this procedure would reduce the operational error to 


2/3 those of the experimental system. Thus we can predict that the 


1g 


standard deviation of the operational line-of-position errors for day- 
time will be 0.32 miles and nighttime, 0.82 miles. [3] 

In the case of the non-correlated lines of position, the root-mean- 
square fix error is usually taken as the square root of the sum of the 
squares of the two line-of-position errors divided by the sine of the 
crossing angle between the two lines of position. If the two position 
lines have the same accuracy we can compute a daytime rms fix error of 


0.48 miles and a nighttime rms fix error of 1.22 miles. [3] 
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Section Three 
Determination of Relative Distance Error Distribution 

3-1 Omega Implementation 

At the initial planning stages of an ASW exercise, recorders capable 
of recording Omega information from the receivers which in turn receive 
signals simultaneously from two different pairs of pre-selected Omega 
transmitters will be installed. The tapes of raw recorded Omega data for 
each individual vehicle will be transferred to a computer facility when 
convenient or at the termination of the ASW exercise. The tapes will then 
be processed in order that the navigation information is in proper form 
for reading into a computer program. A computer program, such as the one 
referenced in Appendix B, will reconstruct the vehicle tracks in such a 
manner as to provide the user with relative distance information of 
selected vehicles (i.e. permutations of aircraft and destroyers with sub- 
marines). With this relative distance information, an analysis regarding 
whether or not any two vehicles were within detection range for a partic- 
ular sensor can be accomplished with the aid of a relative distance error 
distribution. 
3-2 Preliminary Determination of a Relative Distance Error Distribution 

The actual tracks along with the estimated tracks of the vehicles 
must be available in order to determine the relative distance error dis- 
tribution for an area. The actual track of a specific vehicle is defined 
as that track which is not in error. The estimated track is defined as 
that track which incorporates the error inherent to the Omega navigation 
system. With the above two tracks defined, the relative distance error 
at any given time is that value obtained by subtracting the distance be- 


tween the actual positions of the two vehicles from the distance between 
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their estimated positions, as shown in figure l. 

It is immediately apparent that the actual tracks of the various ve- 
hicles are not obtainable because a navigator cannot be sure of his exact 
position. Our solution to this aforementioned problem has been to simu- 
late hypothetical ASW tracks and refer to them as the actual tracks of the 
ASW vehicles. The estimated tracks, on the other hand, corresponding to 
these actual tracks can be determined by incorporating known errors from 
a pair of Omega receivers. These Omega receiver errors can be obtained 
by observing the difference between the known fixed position of the re- 
ceiver's antenna and the position observed from the resulting navigational 
Wee nation being received at a given time, as depicted in figure 1. tt 
follows that these actual tracks will be constructed in the general area, 
a reasonable distance from the antenna, and in a corresponding time frame 
with the antenna. At this point it is necessary for our determination of 
the relative distance error distribution to assume that the dominate error 
inherent to the Omega navigational system is in fact the variation in the 
velocity of the electromagnetic waves in the media and is not appreciably 
affected by the mode of travel of the receivers. This assumption is based 
on the fact that the Omega system operates on the phase difference con- 
cept which is not disturbed by the mode of travel, although the signal-to- 
noise ratio may be affected. 

3-3 Simulation inputs 

The simulated estimated tracks make use of Omega strip charts re- 
ceived from the Naval Electronics Laboratories, San Diego, during the per- 
iod of 1-31 March 1963 (see figure 2). Strip charts were used from three 
different receiving antennas located in the San Diego area and two more 


receiving antennas located in the Pt. Mugu and Pt. Huenume areas. This 


18 





ba 


mop 
= I 
Oo 


TAD 
Sas 


‘ 


bee | amb. 


_ 1 mes 


1 
= = a ie 
= —- te 


A 


ee 
- 
‘ 


_-_ 


ee tes 


gives a capability of getting simultaneous Omega navigational information 
from the five different sources to that simulations of various vehicles 
can be achieved. 

To determine the Relative Distance Error Distribution, a computer 
program was written to gather this necessary information using as its 
primary inputs raw Omega data, time and Longitude-Latitude data of the 
hypothetical tracks. It was necessary to physically interpolate the raw 
Omega data from the strip charts. It may be pointed out here that this 
interpolation technique may introduce errors independent of the Omega nav- 
igation system error itself, but the maximum error was 100 yards in our 
interpolation procedure. It can be assumed that if our Relative Distance 
Errors are acceptable, when a more sophisticated technique is employed, 
better Relative Distance Errors should result because of the exclusion of 
these independent errors. 

3-4 Determination of Antenna Errors 

Omega strip chart information for two fixed antenna positions as a 
function of time is used as input to the computer program. The minimum 
interval of time used for the Omega readings was one minute because of the 
"seaman's eye" technique used for acquiring this information from the strip 
charts. Diurnal variation corrections were then applied to the raw Omega 
data. This was accomplished by reading into the program a matrix correc- 
tion table which corresponded to the applicable period of time and loca- 
tion (i.e. month, day, hour, and longitude-latitude). Since this matrix 
varied not only in time but also in longitude-latitude, it was necessary 
to employ a linear interpolation technique to determine the appropriate 
correction for each input reading (see figure 3). As each reading was 


corrected it then went through a longitude-latitude conversion routine 
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SKYWAVE LINEAR INTERPOLATION CORRECTION ROUTINE 
Selected Latitude-Longitude Outlay 
50N-40W 5SON- 36W 50N~32W 

8--(Vehicle's last position) 


46N-40W 46N-36W 46N-32W 


42N-40W 42N-36W 42N-32W 


Omega station pair correction values for Latitude-Longitude 
which enclose Vehicle's position for specified time 


Time-0900 


50 51 


| 


Ly ng nnn), 5 


@ --(Vehicle's last 
Lg position) 50 


Linear interpolation between values 


Time Interpolation 


Time 0900-48 
Specified Latitude~Longitude value 
Time 1000-53 
Time of position desired - 0930 


Time(0900-09 30) X(55-48) Value = 3.5 


POSITION CORRECTION ~- 49.5 + 3.5 = 53 
Figure 3 


Skywave Linear Interpolation Correction 
Rout ine 
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which is outlined in detail in Appendix C-III. Finally the vector error 
was found by comparing this converted reading which represents the estima- 
ted position of the antenna with the actual antenna position so that the 
inherent error associated for that particular antenna can be determined. 
This process was repeated for the entire input readings for both antennas. 

At this point it is possible to determine the relative error necessary 
to determine the Relative Distance Error Distribution of the reconstructed 
tracks if one is completely convinced that maneuvering of the ASW vehicles 
have no effects on the relative error and if possible smoothing techniques 
were not to be investigated. This relative error can be simply found by 
algebraically subtracting the relative distance of the two actual positions 
of the antennas from the relative distance of their estimated positions at 
a given time (see figure 1). | 
3-5 Hypothetical Track Simulation 

To add intuitive appeal to the problem, actual tracks simulating 
typical ASW tracks for earerare destroyers, or sauiseines were developed. 
This also allows the opportunity to investigate the effects of smoothing 
techniques on the overall relative distance error distribution. An aux- 
lliary program (i.e. see Appendix A) was written to accept Longitude-Lat- 
itude information for constructing an “actual track". To these longitudes 
and latitudes at a given time were added the components of the vector 
error from one of the antenna's to obtain the estimated tracks (i.e. for 
example, antenna A's errors construct a submarine's estimated track and 
antenna B's oe construct a destroyer's estimated track as shown in 
figure 1.) 


Smoothing techniques were then employed in an effort to smooth the 


derived estimated tracks. First by the use of the technique of fitting 
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the estimated track by the method of Least Squares Curve fitting with or- 
thogonal polynomials (refer to Appendix C-II), and secondly by smoothing 

smaller segments of the estimated tracks known as the technique of Expo- 

nential smoothing (refer to Appendix C-I). A new relative Distance Error 
was determined by comparing this smoothed estimated track with the hypo- 

thetical track as a function of time. 

The overall computer program not only has the capability of printing 
out all of the pertinent information such as the Relative Distance Errors, 
time, longitude and latitude as shown in Appendix A, but also prints a 
graphical representation of the actual, estimated and smoothed-estimated 
tracks as shown in Figure 4. 

3-6 Results 

Sample points used for determining the Relative Distance Error Dis- 
tribution were generated from the raw Omega data. Partial sampling of the 
data was accomplished by using permutations of date and time intervals from 
the available data. As an example, a day's transition periods (i.e. sun- 
rise and sunset) were included as well as continuous sections of day and 
nighttime readings. The antennas (i.e. those antennas located in San Diego, 
Pt. Mugu, Pt. Huenume area), which provided the raw Omega data, were loca- 
ted so that trials could be conducted on pairs of antennas which were not 
only in close proximity but also of considerable distance from one another 
(i.e. Pt. Huenume and San Diego). 

A cumulative Distribution of the Relative Distance error (defined to 
be the probability that an observed value of the Relative Distance Error 
is less than or equal to some value) was found from using the data which 
included all the trials. The distribution is found to be Normal with a 


mean of ~120 yards as shown in Figure 5. This mean implies that the 
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relative distance for the estimated tracks was on the average 120 yards 
longer than the relative distance of the actual tracks. 

Investigated next was the cumulative Relative Distance Error Distri- 
bution for the relative error sample points associated with those antennas 
located only in the San Diego area. The analysis of this distribution 
showed that normality could not be assumed, as shown in Figure 6, and the 
mean was 3.2 yards. The above two distributions (i.e. Figures 5 and 6) 
cover the minimum and maximum distance for any two vehicles. 

When the smoothing techniques were used on the estimated tracks, the 
Relative Distance Errors were about ten times as large as those found by 
using no smoothing. These results were attributed to the high degree of 
time dependency of the errors associated with the estimated track and the 
inability of the smoothing techniques to follow precisely the vehicle's 
maneuvers. The comparative results obtained are shown in Figure 7. It 
is recommended that no smoothing features be incorporated when two ve- 
hicles are obtaining simultaneous readings from their respective Omega 


receivers. 
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Section Four 
Other Omega Applications 

4-) Introduction 

If the Omega navigational system were employed as the primary data 
source for the reconstruction of ASW vehicles’ tracks, one might expect 
instances in which continuous track information could not be obtained. 
This could be due to mechanical difficulties with the equipment, or, in 
the case of submarines engaged in evasive maneuvers, running at depths be- 
low the maximum depths for Omega reception. This requires investigation of 
techniques, applied to Omega data, to improve absolute navigational accur- 
acies. 
4-2 Omega~NonOmega Relative Error Accuracy 

When only one of the two vehicles can receive Omega information con- 
tinuously, the track of at least this one vehicle should be constructed, 
so that the tracks latitude and longitude are available for short incre- 
ments of time. A representative sample of the absolute errors associated 
with these latitudes and longitudes are shown in Figure 8, under the col- 
umn for non-weighted vehicles. It should be emphasized that these abso- 
lute errors were obtained from data collected in 1963 and present day 
refinements of Omega have predicted that standard deviations of the opera- 
tional line-of-position errors during the daytime will be 0.32 miles and 
nighttime, 0.82 miles. (3] Trials simulating this situation were conducted 
using vector errors to determine the estimated track, which may be depend- 
ent with respect to other errors for that specific track, but independent 
of the errors associated with a second estimated track. One of the esti- 
mated tracks was assumed to be determined by the use of Omega and the other 


estimated track by some other navigational source such as satellite, celes-~ 
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Onera-NonOmega Weighing Technique 
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trial, loran or dead reckoning. Omega-NonOmega trials were used to deter- 
mine if smoothing the Omega estimated ieee ea improve the relative dis- 
tance errors. This would be expected if the errors in the estimated 
position for each vehicle are independent. Figure 9 shows the results 
from such trials and as expected it appeared that smoothing was in fact 

an improvement. It is suggested that the exponential smoothing method be 
used under Omega-NonOmega circumstances because of the fact that there are 
no navigational error correlations between the vehicles. Exponential 
smoothing was favored over least squares because of its ability to follow 
the vehicles maneuvers was more precise when a high smoothing constant, 
say & = .04, was used. In actual practice the second track would be re- 
constructed in such a manner as to obtain the best absolute accuracy pos- 
sible. This may be accomplished by various methods such as the procedures 
now being used in the formal reconstruction sessions or os recording 
continuous navigational information from a secondary source, e.g. ships 


inertial navigation system. 


4-3 Improvement of Absolute Accuracies by Determining an Expected Esti-~ 
mated Position. 


Improvement in the relative distance error can be realized under 

these Omega-NonOmega circumstances if one could reduce the absolute error 
of the Omega navigational ay8eim: Appropriate weighting of other naviga- 
tional information is proposed. An expected estimated position is found 

by using other navigational information to a particular ship in conjunction 
with Omega. This may be accomplished by weighting the obtained latitudes 
and longitudes by some factor, say Aj» where i denotes a particular type 

of navigational system, and the sum of the A'S must equal one. The value 


of A. is obtained by a comparative judgement of the respective accuracies 
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2 


: longitude 
of the navigational sources (i.e. .3 ( ni a ) of Omega plus 


Pe | i ) determined by seh ctiente’ ). From this expected estimated 
Satellite 
position and the estimated position found when using the Omega navigation- 
al system, a vector difference can be obtained. This vector difference is 
then applied to a number of past and future Omega positions. A sample of 
the results obtained when using this technique is displayed in Figure 8 
with & = 0.4. It may be desirable to incorporate a time factor, say B. 
to the vector difference constant in order to take into account a de- 


creasing effect on the vector error. It is recommended that this weight- 


ing technique be considered when Omega-NonOmega situations arise. 
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Section Five 
Operational Program 

The knowledge gained from Sections Three and Four was incorporated 
into the operational program (i.e. the program used at the reconstruction 
sessions) aS shown in Appendix B. In Figure 10 a flow diagram is refer- 
enced which displays the basic features and the structure of the opera-~ 
tional program. 

The implementation of this operational program during the formal re- 
construction of ASW exercises would aeocide relative distance information 
when Omega is recorded continuously by both vehicles. In conjunction with 
this relative distance information, employment of a previously determined 
Relative Distance Error Distribution (i.e. such as the ones mentioned in 
Section Three) can be used. 

The implementation of this Relative Distance Error Distribution will 
be illustrated through an example shown in Figure 11. It is when the 
vehicles come within detection range of a particular sensor that is of 
primary concern. Therefore of interest are those sample points which lie 
in one tail of the Relative Distance Error Distribution, which reflects 
the situation when the actual relative distance is less than the estimated 
relative distance. For this example the Relative Distance Error Distribu- 
tion found for the antennas located only in the San Diego area will be 
used, see Figure 6. 

The mean of this Relative Distance Error Distribution is so small 
(i.e. 3.2 yards), it is assumed for all practical purposes to be zero. 

As shown in Figure 11 positions (a) and (b) reflect the estimated posi- 
tions of vehicles A and B respectively at time 0900. The Relative Dis- 


tance Error Distribution could be superimposed over one of these positions, 
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Operational Flow Diagram 
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say (b), with the mean situated at (b). The distribution then represents 
the sample points of the various estimated relative distance readings. 

It may be desirable to know what the relative distance between the 
ships must be so that, say ten percent, of the Relative Distance Error 
sample points fall outside of the assumed maximum detection range of a 
particular sensor. For this example, relative distance between the ships 
is found to be - 468 yards from the maximum detection range as shown in 
Figure 6. Therefore as shown in Figure 11 if the maximum detection range 


was 1500 yards, then the estimated relative distance was 1032 yards. 


3/ 


Section Six 
Summary 

It was shown that the position coordinates for a particular vehicle 
can be obtained as a function of time and that a continuous track of the 
vehicles may be graphically displayed. This information requires less 
than five minutes of CDC 1604 computer time for 16 hours of continuous 
interplay. 

A methodology for predetermining a relative distance error distribu- 
tion was demonstrated and an actual distribution for the San Diego area 
was established. This was accomplished by using Omega navigation informa- 
tion collected from receiving antennas whose exact positions are known. 

Attempts to smooth the vehicles estimated tracks did not improve the 
relative distance errors. This apparently resulted from disturbing the 
navigational error correlation between the vehicles. However, when abso- 
lute accuracy of navigational positions was desired (i.e. when only one 
vehicle was receiving Omega information) smoothing the estimated tracks 
was an effective feature. Improvement on the absolute accuracy was also 
observed when an expected estimated position was determined by using other 
navigational sources in conjunction with Omega. 

It is recommended that in view of these desirable characteristics, 
that the United States Navy continue the investigation into the use of 
Omega triangulation information for track reconstruction techniques in 


ASW reconstruction sessions. 
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APPENDIX C-1 
Exponential Smoothing 

Exponential Smoothing is a special kind of moving average that does 
not require keeping a long record of previous positions. Like other 
moving averages, it has a stable response to changes, but the rate of 
response can be adjusted readily. 

Suppose we adopted some such rule as follows: To get a new estimate 
of the average position of any vehicle, add to the previous estimate a 
fraction of the amount by which new positions exceeds that estimate. Note 
that under this rule, positions below the estimate is the Same as saying 
it exceeded the estimate by a negative quantity, and adding a fraction of 
a negative quantity would decrease the estimate. 

The fraction used is called a "smoothing constant" and is represen- 
ted by &. We could abbreviate the rule in the form of an equation: 

New estimate = old estimate + &(new position - old estimate) 


a 


or 


New estimate = & (new position) + (1 - & ) (old estimate). 
This above process could continue indefinitely. A simpler representation 
may be achieved if we call the new position Position, (Pi); the next pre- 


vious position Position (Pj )> the position two sequences ago Position, » 


ik 


and so on, with the position k sequences ago being represented by Posi- 


tion, (PL). Our expression now becomes 
New eStimate = AP, + A(C1l -aH% )P, + Al - A)°P, + (1 - x)°P, 


+. t o& (1 - of SP + (1 - sos (estimate made k 


sequences ago). 
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The new estimate is then a sum of positions in past sequences, each multi- 
plied by a definite weight. If the sum of these weights is equal to one, 
then the estimate can be called an average. 

Therefore the basic rule of exponential smoothing is 


New average = (new position) + (1 - &) (old average). 


How to Select & 
The value chosen for the smoothing constant &% determines how much 
of the past positions has any significant effect on the estimate of the 
average. 

If a small value, say & = 0.01, is chosen for the smoothing constant, 
the response will be slow and gradual, since it is based on the average of 
very many past positions. A high value, say & = 0.5, will cause the es- 
timates to respond quickly, not only to real changes, but also to the 
random fluctuations. 

To evaluate alternative choices of a value for the smoothing constant, 
the accuracy of the positions must be considered along with the effects 


due to Auto-correlation. la] 
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APPENDIX C-IITL 
METHODOLOGY TO BE USED FOR COORDINATE CONVERSION [6] 


The procedure to be — for conversions of raw Omega data to abso- 
lute Latitude/ Longitude (or the inverse) will be based on linear inter- 
‘polation between a set of appropriately selected and stored points. 
The accuracy of the procedure will be insured by the methodology 


used to select the points to be stored. 


The basis for the conversion in either direction will be three matrices 
of stored data. The first two define the crossing points in latitude/ 
longitude of integral (but not necessarily consecutive) lanes of Omega 


phase differences. These points can be defined as: 


L,; = longitude of intersection of Omega lane i with 
Omega lane j 
4;; = latitude of intersection of Omega lane i with 


Omega lane j 


where the first index ''i" refers to the lane number assoc- 
] 
iated with phase difference between one set of two 


Omega signals and the second index ''j'' to another. 


In addition to the position in latitude/longitude, it is necessary to 


store the actual Omega lane numbers of these crossing points. 


wr” L% 


Byrn 
(&{,, = lane number associated with the 'r!! stored cross- 
rs th ; 
ing point for the 'k''” station pair. There will 
normally be only two station pairs being considered 
at one time and it will be assumed that 0,, 8d, ms 


to make searching more convenient. Note that for 
k=l, r=i and for k=2, r5j. 


L. Omega to Latitude/ Longitude 
Reference to figure 1 indicates that one problem to be solved 
is to determine the Latitude/ Lorain (L*/24*) of the Omega 
fix (Q¥, 3) where the subscripts l and 2 define the station 


pair as the "k'' index does in the OQ, matrix. 


By some search procedure it can be determined that Q# lies 

between Qo and Qa nan, and similarly, Qs, n < OF < Qa nt+le 
ae 

In addition 0) is closest to 034 and {3 is closest to yxe The’ 


equations defining L* and t* then become: 


' Co i = lis ) “ 
L¥ = Log 4 et =e ACQ 1 
’ Ae ; i (1) 

L ant La ) nS 

cE oe AQs 
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FIGURE 1: Geometry of conversion problem 


, a ee » (Lpta nn - n) 
4% = 43,9 + Mo anal acy 


Leiner ~ 4 ) te 
v Aa. AQs 


where AQ, = (82a) ata - Ga, 2): 
AQs (Qe nea = Qa, 2) 
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